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(yQ ' Abstract 

cn '. 

\^ . The readout noise of a H2RG HgCdTe NIR detector from Teledyne 

cn I is measured at a temperature T=110K. It is shown that a Fowler 

lO ' mode with n = 240 allows to reach a noise of 2.63e (single read). A 

description of the power spectrum in terms of 3 parameters reproduces 
the variation of the noise as a function the number of Fowler samples, 
as well as its dependence on the periodicity of the sampling. The 
variance of the noise decreases with frequency with an effective power 
of 0.62 in our measurement domain. The behaviour of the detector 



C^ ' under different experimental conditions can then be predicted. 



1 The Apparatus 

The measurements described in this paper were carried out in a dedicated 
setup built to evaluate Hawaii 2RG (HgCdTe) detectors from Teledyne. The 
detector was on loan from LBNL in view of the evaluation of its performance 
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when used in a spectrograph for the JDEM project [1]. The cryostat can 
be operated in a range of temperature extending from llOK to 160K with 
fluctuations smaller than 0.1 K, and its equilibrium temperature in the ab- 
sence of heating is about llOK. The polarisation of the substrate was chosen 
as Vsub — Vreset = 0.4 V in all the data analysed in this work. Additional 
experimental details are provided in a previous paper [2] where the conver- 
sion factor of our setup (e/ADCU) at T = llOK was evaluated to be 2.042 
e/ADU for a single pixel. The goal of the present measurements is to lay the 
ground for a determination of the frequency distribution of the noise, so as 
to be able to predict the behaviour of the detector in different experimental 
conditions (Fowler samplings, interval between groups, etc.). The merit of 
the power spectrum analysis is also the link it provides with the physical 
processes in the CMOS semi-conductor, such as the trapping and detrapping 
of electrons with a distribution of time constants. 

2 The measurement method 

In the 2Kx2K H2RG detector, a window of 31x31 pixels has been selected. 
The clock frequency is lOOkhz, as anticipated in the SNAP- JDEM project, 
and the time needed for the readout of 1 frame of 30x30 pixels is 5 = 17.5ms . 
The non destructive frame readouts are organised into 'groups' of 250 frames, 
read at the clock frequency of lOOkHz and separated by a time interval which 
can be tuned, with the clocking of the pixels stopped. Taking into account 
the readout of the frames, the periodicity of the readout of a given pixel 
between 2 consecutive groups varies from 5s to 13s. The selection of a small 
window allows to reach a higher repetition rate for the frame readout, and to 
save on the overall measurement time as well as on computing resources. All 
the frames are stored on disk for further analysis. The noise is characterised 
by two variance measurements: 

-the frame to frame variance, with the readout interval oi 6 = 17.5ms for a 
'window' frame 31x31 as considered here, is convenient for a calibration of 
the detector and its readout. It was studied in [2]. 

-the Fowler(n) group to group variance measured by the spread of the dif- 
ference of the averages of 2 consecutive groups for a given pixel will allow to 
characterize the noise performance of the H2RG on longer time scales, and 
will be investigated here. 

The Fowler (n) noise is evaluated in each pixel from the differences D^ of 
the average signal in groups k and k — 1. 



-Dfc = - 2^ S{to+kA+iS) - S(to+ik-l)A+iS) (1) 



The number of frames is 250 in the acquisition, but will be varied from 1 
to 240 in the offline analysis as the first 10 frames of each group are always 
ignored. The noise is measured by the average of Z}„(A)^ over the A^ + 1 
groups stored during the exposure. The number of groups in the on-line 
acquisition was actually 200 in the exposures studied here, but the first 10 
groups were eliminated from our investigation. 
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The subtraction of the 
reference channel, which 
follows synchronously the 
level of a fixed capaci- 
tance, and corrects pos- 
sible bias changes reduces 
the variance when the num- 
ber of Fowler samplings 
n < 10, as shown in Fig- 
ure 1 for a group period- 
icity of 5.23s. For larger Figure 1: Noise reduction from reference channel 
numbersof reads, the vari- and common mode subtractions. The second one 
ance of the noise is how- is applied on top of the first one 
ever slightly increased, by 

typically O.lADU^. We attribute this to low frequency voltage offsets, which 
are at least partially cured by the next step. 
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3 Common-mode subtraction at low frequency 

The reference channel subtraction takes out most of the high frequency com- 
mon modes, but we still observe the impact of common modes remaining 
at lower frequencies in Figure 2a). This figure shows the remaining corre- 
lation between the differential flux observed for 2 randomly selected pixels 
between 2 consecutive groups (in the absence of any illumination), and the 
differential flux of the sum over all pixels. Both are averaged over the 240 
Fowler acquisitions of the group, and the 200 groups of the exposure are 
shown. The value of the slope of the correlation is stored for each pixel, 
and its distribution is shown for all pixels in Figure 2b), where it is seen to 
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Figure 2: a) group to group correlation between pixel flux (Fowler=240) 
and full flux: 2 pixels have been selected at random, b) distribution of the 
correlation slope (Fowler = 240)found for all pixels 



be compatible with unity within the measurement errors from the noise. A 
'common' offset for each group is then obtained from a second straight line 
fit throughout the 'full' (integrated over all pixels) observed fluxes of all 200 
groups as a function of the group number. The resulting variance reduction 
is shown in Figure 1 as a function of the Fowler number n, and the resulting 
distribution of the 240 corrected flux differences D^ is shown in fig 3, for a 
typical exposure at 11 OK with Fowler samplings n = 240. The rms of the 
distribution is 4 electrons. As this value arises from the difference between 
the averages over 2 (consecutive) groups of 240 frames, it can then be con- 
sidered that the noise in a single Fowler readout is 2.85 electrons. We shall 
see in section 5 that part of this noise can be assigned to a parasitic photon 
flux, and that the actual performance reached is better. The improvement 

mean of 240 trames group to group differences | 
800 = 




Figure 3: group to group difference distribution for Fowler(240) 



observed with this common-mode subtraction may be a combination of in- 
trinsic detector readout properties and of the specific implementation of the 
outside low voltage supplies. 



4 Evolution with the number of samphngs 

The evolution of the variance of the difference Dk between consecutive groups 
as a function of the Fowler(n) number is shown in Figure 4 for an intergroup 
periodicity of 5.23s after the reference channel and common mode correc- 
tions have been applied. The duration of the full exposure (250 groups of 
250 frames each) is 1280 seconds. The value reached for Fowler(240) corre- 
sponds to a noise of 2.85e for a single group read. It is apparent that the 
variance does A^'OT decrease like 1/n, as would be expected for independent 
measurements,which would lead to a group to group noise of 1.6e at n = 240. 
There are two main contributions to these correlations: the non destructive 
readout itself, and the time domain correlations induced by the trapping and 
detrapping of electrons. This behaviour has been noted in several previous 
publications, such as [3], [4], [5], and implies the presence of time dependent 
correlations, which can be described by a frequency power spectrum. The 
connection between the time constant of trapping sites and the frequency 
power spectrum was first stated by [7], and a recent presentation can be 
found in the lectures of [8]. 

5 The extra flux 

Before turning to the frequency power spectrum, we must take care of the 
other contributions to the observed variance. Even in the absence of any 
illumination inside the cryostat. The ADC level rises linearly with time, at 
the rate of 0.6 e/pix/s whatever the group periodicity. This rate is not a 
leakage current, as it is (almost) insensitive to the operating temperature, but 
it could be due either to a parasitic photon source or to some baseline drift. 
We shall assume here that the cause is a photon source subject to stochastic 
fluctuations, and we consider that the quantitative agreement with the data 
reached under this assumption argues in its favour. The impact of such a 
constant photon source on the Fowler(n) variance is recalled in Appendix A 
for completeness, although it has already been evaluated in the litterature 
in a slightly different context, such as in [6], where the impact on a straight 
line fit is directly evaluated. The result is of the form 

{PD =DI + fdif (4) 

where (P^) is the flux contribution to the variance, DI the intergroup contri- 
bution (the flux already included in the n Fowler frames must be omitted), 
dif the frame to frame contribution, and the factor / takes into account the 
correlations arising from the non destructive readouts, as shown in Appendix 
A. 

^ 2dif n(n -l)(2n-l) ,, 

n-^ 6 



6 Frequency spectrum of the noise 



We show in Appendix B how to relate the Fowler averaging to the auto- 
correlation products. As the autocorrelation of the noise is itself directly 
related to the frequency power spectrum /(w), the group to group noise for 
any number of Fowler samplings can then be expressed in terms of the power 
spectrum: 
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The power spectrum f{uj) will be assumed to be of the form 



f{u) = A + 
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u^ 
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If this functional dependance were exact, the white 'thermal' amplitude A 
would be a positive number, and a might generaly be expected to be close 
to 1. Our measurements are however performed in a frequency range where 
the white noise is subdominant, so that A is just a parameter of the fit, and 
it turns out it will be found to be small and negative. The power a found 
is also an effective power in our frequency range, which may combine the 
detector and FET effects. The merit of this 3 parameter formalism is that it 
describes the shape of the Fowler (n) function with an accuracy of 2% and 
has a predictive power for different group or frame periodicities as well. 

As shown in Appendix B, the measured variance (D„(A)^) can be ex- 
pressed as a sum of autocorrelation products, and an explicit expression 
of the integral relating the power spectrum f{uj) to the measured variance 
(-D„(A)^) is derived in the appendix: 
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In this relation, 6 and A are the frame to frame and group to group 
periodicities. The predicted variance for given values of A,B, and the power 
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a, are evaluated at 10 different Fowler samplings (to spare computer time) 
and they are adjusted to the data for 2 different group periodicities of 5.23s 
and 13.11s by minimizing the ~x^ . Given the very high statistical accuracy 
of the measurements, and their reproducibility, the typical computed errors 
on each data point in Figure 4 is 3 to 4 10~^. On the other hand, the 3 
parameter model used for the power spectrum reproduces the data to an 
accuracy of 1 to 2% for the best fits. The errors were rescaled to match this 
difference and follow the observed trend, as seen in Table 1. 

The best fits differ for the 2 data samples, although the power spec- 
tra found are close to each other. The estimated increase in the noise at 
Fowler=240 is shared equally (in terms of variance) between the stochastic 
contribution and the 1/f contributions. 

The (dominant) B contribution amounts to 27.107e/(Hz)^^'^ (5s inter- 
group) and 29.096e/(Hz)~^/^ (13s intergroup). The result of the fit is seen 
in Figure 4a, and the estimated errors on A and B shown as contour plots 
in Fig 4b) are only indicative, as the errors are just set so that the y^ jNdf 
would be of the order of unity. 

Pb vs A, 1 and 3 a contours, a=0.62 | 
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Figure 4: a) Fowler (n) noise measurement and adjusted frequency spectrum 
at a group periodicity of 5.23s b) contour plot of the A and B errors 



7 Prediction of the group periodicity depen- 
dence 



We can also check how accurately the evolution with the group period is pre- 
dicted from the restricted lever arm in frequency provided by the Fowler(n) 
dependence of the noise at a single delay. We show in table 3 the depen- 
dence of the variance expected with 240 Fowler samples on the intergroup 



Table 1: Variance evolution with the number of Fowler samplings) 





Intergroup 


periodicity 5.23s, 


exposure time = 1280s 




Fowler nb. 


Adjusted noise 


Meas. noise 


Meas. error 


Assigned error 


(n) 


(e) 




(e) 




(e) 


(e) 


2 


18.58 




18.28 




0.135 


0.37 


4 


14.37 




14.61 




0.055 


0.29 


5 


13.29 




13.58 




0.033 


0.27 


10 


10.57 




10.81 




0.022 


0.22 


20 


8.55 




8.62 




0.025 


0.086 


50 


6.56 




6.51 




0.045 


0.065 


100 


5.37 
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4.52 
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4.04 
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periodicity 13.11s, 


exposure time = 3275; 
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Fowler nb. 


Adjusted noise 


Meas. noise 


Meas. error 


Assigned error 


(n) 


(e) 




(e) 




(e) 


(e) 


2 


19.04 




18.72 




0.025 


0.37 


4 


15.02 




15.26 




0.032 


0.30 
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13.98 




14.27 
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0.28 


10 


11.35 




11.56 
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20 


9.40 
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50 


7.43 




7.35 
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6.26 
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150 


5.74 
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5.54 




5.55 
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200 


5.37 
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240 


5.15 
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0.052 



delay. The stochastic Poisson contribution from the assumed 'extra flux' is 
substantial, and the agreement with the observations suggests that the gen- 
eral features of the physical contributions to the noise are understood. A 
small discrepancy is observed at a group periodicity of 13s, confirmed by an 
inspection of the full Fowler prediction in Figure 4a). We are tempted to 
assign this effect to a frequency dependence of the power a, which will be 



Table 2: Adjusted power spectrum parameters at 5.23 and 13.11s group peri- 
odicities 
Exposure time/ Thermal 1/f 1/f power x^/Ndf 



period(s/s) A (/iV^Hz-^) B (/xV^Hz 



2U^-1^ 



a 



1280/5.23 
3280/13.11 



8.445 
-1.525 



734.819 
846.563 



0.62 

0.57 



6.063/^ 
6.000/^ 



Noise Spectrum 




Figure 5: Frequency power spectrum of the noise with the adjusted param- 
eters at group periodicities of 5.23 (continuous line) and 13.11s (dashed) 



investigated in a later work. 



Table 3: Variance evolution with group periodicity for Fowler(240) 
Period/Exp. time Var. (e^) Var. (e^) Var. (e^) Var. (e^) 
(s)/(s) (flux) (frequency) (predicted) measured 



5.23/1280 


2.396 


13.764 


16.160 


16.646 


6.578/1315 


3.427 


14.949 


18.376 


18.768 


7.304/1461 


3.941 


15.452 


19.393 


19.564 


8.756/1751 


4.905 


16.085 


20.991 


21.328 


10.208/2042 


5.956 


16.667 


22.623 


23.088 


13.112/3275 


8.066 


17.472 


25.538 


27.041 



It is seen that the actual readout contribution to the group to group vari- 
ance is 13.8 e^, i.e. 3.71 electron for the difference between 2 Fowler(240) 
readouts and a group periodicity of 5.23s . The extrapolation of the pa- 
rameters adjusted at an intergroup period of 5.23s is compared to the full 
Fowler(n) curve at a period of 13.11s in Figure 4a. It is seen that the pre- 
dicted noise slightly undersestimates the measured one, as also observed in 
table 3. The adjustment at a period of 13.11s has an unphysical negative 
value of the white noise component, as seen in Table 2 and Figure 5a). This 
feature has a minor impact as all the experimental data is carried in a fre- 
quency range smaller than 370 Hz, where the 2 frequency spectra are close 
to each other. 

8 Conclusion 

We have shown that the combined use of the reference channel and a common- 
mode correction lead to a significant reduction of the detector noise. The 
second correction can be performed with the help of non illuminated pixels 
(unavailable in the window mode) in the standard operation of the H2RG 
detector. With the H2RG lot # 40 which we measured, we reach a noise of 
2.63e (single read) with Fowler(240) in a 1250s exposure. This is degraded 
to a noise of 2.96 e for a 3280s exposure with an intergroup periodicity 
of 13s. A description of the power spectrum of the noise in terms of 3 
parameters also accounts quantitatively for the Fowler(n) variation of the 
noise at different group periodicities within an accuracy of 2%. The power 
law exponent found to reproduce the data is about 0.62. This 'effective' 
power law is only approximate, and the impact of additional parameters 
improving the description of the power spectrum will be investigated in a 
later work, as well as the effect of different polarisation voltages. 
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Appendix 
A Flux and leakage current contribution 

S and A will be the interval betweeen frames and between groups. In the 
Fowler mode, an average is taken between signals at t,t + S,t + 2S, ...t + 
{n — 1)5 and the difference is formed with the homologous signal values at 
t + A, t + A + 5, t + A + 25, ...tA + {n — l)5. The contribution of the stochastic 
fluctuation is evaluated next DI and di will be the mean intergroup and 
frame to frame fluxes, S{DI) will stand for the intergroup signal fluctuation, 
and 6i{k) will be the flux fluctuation of group i at frame k (the fluctuation 
of the difference between frames k and k — 1). 



S^ = DI + 6{DI) 

52 = DI + 6{DI) + di + 62(2) - 6i{2) 

53 = DI + 6{DI) + 2di + 62{2) + 62{3)-6i{2)-6i{3) 

Si = DI + 6{DI) + 3di + 62(2) + 62(3) + 62(4) -Si{2)-6i{3)-6i{4) 

Sn = DI + 6{DI) + {n-l)di + 62{2) + ...62{n)-Si{2)-..Mn) 

The average is formed 

Oi + O2 + 03 + ...On 



s 

n 
The fluctuating Poisson contribution to the result is 

-((52(2)-5i(2))*(n-l) + (52(3)-5i(3))*(n-2) + (52(4)-5i(4))*(n-3)+...+((52(n)-5i(n))*l 
n 

The variance of 52{i) and 5i{i) are all equal to the number of electrons 
between 2 frames dif (in number of electrons) 



1 
n 
where the sum of the n — 1 first squares appears 



5{Sf = DI+ ^2dt}{l + 2^ + 32 + ....(n - 1)2) 



,2 , o2 , q2 , ,^ 1^2 _ (^- !)("-) (2^-1) 



1^ + 2^ + 3^ + ... n-1 



6 

The stochastic contribution to the variance Fowler(n) variance will be (in 
electron units) 

.(5f^o/ + l2..;'»-^><"'P»-^> (9) 
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To be added quadratically to the frequency spectrum noise. It is seen that 
the stochastic error on the group average grows hke the number of Fowler 
averages, and will exceed the readout error at some point. High Fowler 
numbers harm at large fluxes. 

B The Wiener relation 

Time domain correlations of the noise are induced by the interplay of the 
physical properties of semi-conductors. As shown by Wiener, they can be 
described by the frequency power spectrum of the noise,which happens to 
have a quasi universal behaviour as a function of frequency, of the form 
A + B/u°'. Assume 

x{t) = / doo f {uj)e—iujt 



1 r+T 



<x{t)x{t + A)> = lim — / dtduJiduJ2e~''^^'e''^^'+^f{uJi)f*{uj2) 

T^-oo 21 J —T 



hm — / rfa;e-^/(c^)/*(a;) 



As the result should be real 

f{-uj) = /M 



< xit)xit + A) >= -( / dcufiu)fiure'-^ + / dufioj)fiuye^-^) 

^ Jo J-oo 

So that taking advantage of the reality relation 

POO 

< x{t)x{t + A) > = - / rfw(e^^^ + e-'^^)f{uj)f{ujy 

T Jo 

= — y ducos{coA)f{u)f{ur 
including the factor 271 /T in the normalisation of the power spectrum / 

< x{t)x{t + A) >= dco{cu)cos{cuA)f{uj)f{ujy (10) 



JO 

The shape of the power spectrum | f{oo) p will be assumed to be of the form 

where A is the intensity of the white noise component and B related to the 
strength of the l/u (1/f) contributions. 
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C The group to group averages 

The measured quantity is the variance group to group differences in the 
Fowler(n) averages. We introduce the time intervals S between 2 reads of the 
same pixel in consecutive frames, and the equivalent interval A between two 
readouts of the same pixel in the same frame in consecutive groups. 
For the consecutive groups k and k + 1 

Dk = -ix{tk+A)-x{tk)+x{tk+5+A)-x{tk+5)+x{tk+25+A)-x{tk+25)+etc...) 
n 

(11) 
<Dk>=^Y.^^ (12) 

k 

is the time average of D (i.e. the fluence). The noise variance will be defined 

by 

<Dl>=^Y.^D,-<D,>f (13) 

fc 



<Dl>= 

— (2n < x^ > -2n < x{t)x{t + A) > 

+2 < x{t + A)x{t + A + 5)> +2<x{t + A)x{t + A + 2(5) > +... 

-2 < x{t + A)x{t + 6) > -2< x{t + A)x{t + 25) -2 < x{t + A)x{t + 35) > -.. 

-2 < x{t)x{t + A + 5)> -2< x{t)x{t + A + 25) <-2< x{t)x{t + A + 35) > - 

+2 < x{t)x(t + 5)> +2< x(t)x{t + 25) > +2 < x{t)x(t + 35) > +...) 

(n — 1) similar terms 

+2<x{t + A + 5)x{t + A + 25)> +2<x{t + A + 5)x{t + A + 35) > +... 

-2 < a;(t + A + 5)x{t + 25) > -2 < a;(t + A + 5)x{t + 35) 

-2 < a;(t + A + 5)x{t + 45) > -... 
-2 < x{t + 5)x{t + A + 25) > -2 < x{t + 5)x{t + A + 35) > 

-2 < x{t + 5)x{t + A + 45) > -... 
+2 < x{t + 5)x{t + 25) > +2 < x{t + 5)x{t + 35) > +2 < x(t + 5)x{t + 45) > +. 
{n — 2) similar terms 
{...)The summation is continued in the same way 



The correlated products can be explicitly expressed in terms of the power 
spectrum as in the previous section 
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<x^ > 



JMA + ^J 



<x(t)x(t + A)> = I cos{ujA){A + —^ 



And using the time translation invariance 



cos{uj5){A+ — ) 
cos(a;(A-5))(A + — ) 



etc.... Bringing together the correlated products with the same value 

<Dl> = ^jduj{l-cos{uA){A+^) 

+ {n — 1){2cos{uj5) — cosa;(A — 5) — cosa;(A + 5)) 
+ {n- 2){2cosu26 - cosw(A - 26) - cosa;(A + 26)) 
+ (n — 3)(2coscj35 — coscj(A — 35) — cosu{/S. + 35) 
(...) the summation is continued in the same way 

So that 

<Dl>= 

^ r R 

— [2n / duj{l - cosuj/^){A + -^) 

/JD 
duj2cos{uj6){l - cosujA){A + —^) 

du2cos{2u6){l - cosu/\){A + -^) 

+2(n - 3)( / duj2cos{?,u6){l - cosujA){A + — ) 



(14) 

It is straightforward to provide an analytical expression of the sums in 
equation (14). 
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D The summed results 

The summation in equation (14) involves to different series. 

cos{uj5) + cos{2uj5) + ... + cos{{n — \)uj5) 



= R-{^'''' \l,..s } (15) 

Similarly 

cos{ijj6) + 2cos{2uj6) + ... + (n — l)cos{{n — l)u6) 

1 d 

= -——(sin(ijj6) + sin(2u6) + ... + sinin — 1)oj6) 

oco 

1 (9 
6 doo 






(16) 



The Fowler(n) averages expressions can then be derived 

E Explicit result for the real part 

(1 - e"^^) = I- cosuj5 - isinwS = 2 sin^{u6/2) - 2isinu6/2cosu6/2) 
(1 - e*^^) = -2ism(a;(5/2)e^^^/2 



1 _ (,iu>& -2isin{uj5/2)e'^^/'^ 2 sin{u6/2) 
The cosine summation is 

n sin[ujo/2) 

As 

(1 - e'"("~i)^) = -2isin{uj{n - 1)5 /2)e'''^''-^^^'^ 

The real part contribution to the variance < D\ > is 

i_cos(n.S/2)srn((n - 1)^^/2) _ B_ 

n sin{oo5/2) ^ '^ u'^' ^ ' 
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F Explicit results for the imaginary part 

n^6du 1 - e^'^^ 

Aid -2ism(w(n - l)(5/2)e^'^("-^)''/2gi^5/2 

n? 5 doo —2isin{uS/2) 

4 1 d sin{uj{n ~ 1)5/2 ^^■^s/2^ 
n? 5 diij sinoo5/2 

This contribution to < D? > is 



A I d sin{u{n - 1)5 2) . , ^.,,, ^,,^ B. 

n^5du sinu5/2 ^ ' "^ '^ w"^ 

After evaluating the derivative, the imaginary part contribution is 

9 R 

\sir^{nuj5 /2)-nsin{oj{2n-l)5 /2)sin{uj5 /2)){1-cosuj/\){Aa ) 



71^ sin'^{uj 5/ 2) id' 

(18) 

G The final formula 

Adding the real and imaginary contributions obtained in eq. (17) and (18) 
<Dl> = fdu{l-cosuA){A+^)[- (19) 



4 cos{nuj5/2)sin{{n — 1)uj5/2) 

n sin{u5/2) 

2 

+ -^ o/ r/ A sin^(nu5/2) - nsin(u(2n - l)5/2)sin(u5/2))] 

n^sin^[uj5/2) 



17 



B vs A, 1 and 2 a contours, a=0.62 



N 



37.5- 






36.5- 




35.5 



A (nVA/Hz) 



